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We consider issues related to the conformal mapping between the Einstein and Jordan frames 
in f(R) cosmology. We consider the impact of the conformal transformation on the gauge of a 
perturbed system and show that unless the system is written in a restricted set of gauges the 
mapping could produce an inconsistent result in the target frame. Newtonian gauge lies within the 
restricted group but synchronous gauge does not. If this is not treated carefully it could in principle 
contaminate numerical calculations. 



CN " I. INTRODUCTION 

O- 

. Extended gravity theories, where the Einstein-Hilbert Lagrangian density £eh = R is replaced by a more general 
function including terms of higher-order in derivatives of the metric (i? 2 , R^R^, R a ^iivR OL ^ llu • • ■) and couplings to 

O . new dynamical degrees of freedom, have long been of interest in relativity. The last few decades have seen increasing 
<^ applications of these models to cosmology. See [l|, Q and their references for an overview and further details on these 
models and their applications to cosmology. 1 In the last decade attention has focused on exploiting extended gravity 
to model dark energy without the need to introduce exotic particle species. A model frequently employed in this 
context is f(R) gravity (see for example [1, H| for recent reviews), where the Einstein-Hilbert Lagrangian density is 
replaced with an arbitrary function of the Ricci scalar, C = f(R), while the matter couples minimally to the metric 
and follows its geodesies in free motion. This representation of an extended gravity model is known as the "Jordan 

$H . frame" . 

The action can be transformed to a variety of different forms. One of the most common transformations is into 
| the so-called "Einstein frame" where the action is manipulated to isolate a Ricci scalar of a new metric. The residual 
£Xj terms can be interpreted as an effective scalar field to which matter is non-minimally coupled and deflected from the 
£> 1 geodesies of the metric. The field equations are otherwise those of standard general relativity. We review the model 
I/"*) . in the Jordan and Einstein frames in SjTTJ 

Transforming from the Jordan frame to the Einstein frame is extremely useful in the study of f(R) gravity. We 
employ the metric formulation, in which the action is varied with respect to the metric alone, and in which the 
field equations are fourth-order. In the Einstein frame, as in standard GR, the theory is second-order - a significant 
simplification. Aspects of the transformation have been controversial for some time (see for example an d further 

references in . The discussion has centred upon the nature of the equivalence, and authors can be separated [l|, || 
into two camps: those who feel the equivalence is "physical" and observables can be calculated in either frame, 
• • ■ and those who feel the equivalence is mathematical in nature and that observables should be calculated in a chosen 
"physical frame" . We briefly discuss this issue in 811111 

Modern cosmology is the study of Friedmann-Lemaitre-Robertson- Walker (FLRW) spacetimes, perturbed to linear 
or higher orders, and this provides an illustration of the major motivation for our paper. Numerical studies targeting 
the cosmic microwave background (CMB) and the matter power spectrum (such as [9Hl3|) often employ modifications 
to Boltzmann codes such as CAMB [lij]. While it is common to employ a general parameterisation [Tol - [T3| in the 
Jordan frame, working in the Einstein frame (e.g. [j| [r|) provides a flexible alternative. Exploiting the Einstein 
frame allows us to consider any modified gravity model which possesses an Einstein frame and for which background 
solutions in the Jordan frame might be extremely difficult to find, either analytically or numerically. A recent example 
is the study of the string gas cosmology (see e.g. [Hj]), in which the primordial power spectrum is frequently calculated 
in the Einstein frame as opposed to the string frame, as in [T3, IH| • The use of the Einstein frame allows us to import 
our intuition and understanding from standard gravity - at least while we remain within the frame - and, perhaps 
more importantly, it allows us to directly leverage well-tested codes developed for general relativity and which contain 
a wealth of physics and arbitrary collections of fluids. Codes of comparable breadth would be extremely lengthy to 
implement in the Jordan frame, and their fourth-order nature has consequences for speed and stability. 
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1 Recently [|| presented the most general action guaranteed to be free of ghosts when quantising perturbations on a classical vacuum; a 
viable model of extended gravity will be constrained to conform with this action. 
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Considering the nature of the transformation for a perturbed model reveals an issue which to the best of our 
knowledge is as yet unappreciated. Perturbed systems in relativity exhibit the gauge issue - the mapping between the 
fictitious background spacetime and the physical perturbed spacetime is non-unique. In cosmology, gauge freedoms 
allow one to eliminate four perturbative degrees of freedom (two scalar and two vector, where the classification is 
based on how the modes transform on a 3-surface), which specifies the slicing and threading we have chosen for our 
foliation. The transformation between the Einstein and the Jordan frames tangles this choice. We consider the impact 
of this and how one can resolve the resulting gauge ambiguities. We illustrate with a simple f(R) model in a vacuum 
FLRW spacetime, but it should be emphasised that this issue will in principle occur in any perturbed spacetime and 
a wide range of extended theories of gravity. 

We choose to work in the standard ("metric") approach to cosmologic al p erturbations (e.g. [Ill [2(|) and leave 
potential issues in the 3+1 gauge-invariant and covariant (GIC) approach |2l| - [23l | for future study. 2 While many au- 
thors working on perturbed spacetimes in modified gravity choose to employ "gauge-invariant" variables (for example, 
[l9l [25j ) or work in the Jordan frame (such as [26l - [33 |). the analysis of gauge problems is necessary for three reasons. 
Firstly, fixing a gauge (as in [tj 13314351 ] ) . is a valid approach, and any problems when transforming to the Einstein 
frame must be understood. Secondly, gauge-invariant variables are themselves based on variables in a particular gauge 
[36( | . The frequently-employed Bardeen potentials are, for example, gauge- invariant generalisations of the potentials 
in Newtonian gauge. If there is an issue with the choice of the "natural" gauge of a gauge-invariant perturbation, the 
perturbation itself cannot necessarily be trusted. Thirdly, our study is motivated not least by the use of Boltzmann 
codes in the Einstein frame. By far the widest used Boltzmann code is CAMB [14j, which is written in a frame 
equivalent to synchronous gauge. 3 If there is an issue with the frame transformation and synchronous gauge, then the 
initial conditions for modifications of CAMB could be suspect, as could the final results even after being transformed 
back into the Jordan frame. Worse, if there is an issue the error is compounded: each transformation between frames 
will induce errors. Given that the use of CAMB (or similar codes) is of interest, the study of transformations of 
systems in synchronous gauge is forced upon us, despite its known shortcomings: synchronous gauge is not fully fixed 
[3a . l42j and ambiguities in the choice of threading remain. 

The alternative formulation of f(R) gravity employs the Palatini approach, where the metric and the connections 
are treated as independent objects. The Palatini approach is popular in cosmology (e.g. [H, 0, Hi| and their references) 
and often employed in studies of cosmological perturbations (e.g. (45l - l47| ). The nature of the transformation between 
the frames implies that the issues we highlight in this paper are also relevant in the Palatini approach. Since the field 
equations in the Palatini approach are intrinsically second-order the usefulness of the transformation is somewhat 
lessened, but if a system is more easily evaluated in the Einstein frame then the potential ambiguities should still be 
considered. 

The gauge problem is discussed in ^IVI A and illustrated in £|IVB . We close the paper in <JV]with a brief discussion. 

We employ a signature ( — h ++) and the Ricci tensor R^ v = R a miJ = r" a + . . . (as in for example d, HH, l49j ) 
and generally follow the notation of Q. An overdot will represent differentiation with respect to conformal time 
dr) = dt/a, and the conformal Hubble rate is written T-L = a/a. 

II. F(R) GRAVITY IN THE EINSTEIN AND JORDAN FRAMES 

A. f(R) Gravity in the Jordan Frame 

In f(R) theory the Einstein-Hilbert Lagrangian density is replaced with a general function of the Ricci scalar, leading 
to the action 

S = ^ fd 4 xV=gf(R) + [d 4 xC M (g MV , *m) (l) 



2 The issue of perturbed models of modified gravity in the 3+1 GIC approach has been considered in [24H . The authors wish to thank 
Sergei Odintsov for bringing this work to our attention. 

3 Of the other modern Boltzmann codes, cmbeasy [37ll can be used in synchronous or Newtonian gauges, and CLASS 38] is best tested in 
synchronous gauge although it was updated with support for Newtonian gauge after the release of the first version of this paper. Very 
recently, CosmoLib 13911 was released, programmed entirely in Newtonian gauge. COMICS 14011 could be employed in both Newtonian 
and synchronous gauges, but was rapidly superseded by CMBFast 41] which could only be used in synchronous gauge. Both of these 
codes are now obsolete, with CMBFast in turn superseded by CAMB and cmbeasy. Our comments concerning CAMB and synchronous 
gauge equally apply to other Boltzmann codes in other gauges. 
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where denotes the matter fields. The field equations follow by varying this action with respect to the metric, 

FR, V - \f9»» - V,V,F + g^V a F = n 2 T^\ = ~ J^f g » v - (2) 

Here F = df /dR. The matter stress-energy tensor obeys the standard conservation law, 

VI^ = 0. (3) 

This form of the theory is known as the "Jordan Frame" . By introducing an auxiliary field it can be shown (e.g. 0,111) 
to be equivalent to a Brans-Dicke theory with wbd = so long as d 2 f/dR 2 ^ 0. In the Brans-Dicke representation 
the theory has also been described and studied as an example of "extended quintessence" (e.g. [5l[ and its references). 
The flat FLRW metric in the Jordan frame can be written 

ds 2 = a 2 (r/) (-dr? + 5 lJ dx i dx J ) . (4) 

For more details on f(R) theory, see for example 0, H3, [Elf and their references. 



B. f(R-) Gravity in the Einstein Frame 

f(R) gravity can be cast in a more familiar form through a conformal transformation 

9v» -> 9^ = ti 2 g^ gT = ^V, g = Q 8 g. 



(5) 



A conformal transformation leaves null geodesies invariant, implying that the causal structure of spacetime is unaf- 
fected by this mapping. Writing In Q — oj, the Christoffel symbols, covariant derivatives of a scalar and a covector, 
the Ricci tensor, Ricci scalar, stress-energy tensor, stress-energy conservation transform [l|, 0, 0, as 



R^u = R„u - 2V M V„w - g^V a V a u; + 2V (U wV I/ u; - 2g^V a ioV a to, R = fl~ 2 (R - 6V^V M u; - eVwV^w) , (6) 

rjn{M) _ _2 T (M) 



and the matter is deflected from the new geodesies: 



V M- t (M) 



-7 1 (M)V !/ W. 



Writing / = RF + (/ — RF) = RF — U, the action |T]) can then be rewritten as 

F 



S 



2n 2 n 2 



R + 6V^V M w - 6VwV M w ) - 



u 



Cm {a- z g„ v ,^ M )^x. 



(7) 



(8) 



The Einstein frame is then defined by isolating a Ricci scalar in the action to serve as the Einstein-Hilbert Lagrangian 
density; this can be found by choosing 



Q = F 



(9) 



which requires the condition F > to be always satisfied. Casting the rest of the action into the form of a scalar field 
minimally-coupled to gravity, 



produces 



S 



(10) 



(11) 
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This action describes general relativity in the presence of a scalar field minimally-coupled to gravity but non-minimally 
coupled to matter - an example of a "coupled quintessence" (e.g. [5ll l54U58j | ) model. 4 The stress-energy tensor and 
potential of the scalar field are 

T$> = v>V„«£ - ig^V^V^ - V($)g^, V{4>) = ^f^- (12) 
and the field obeys the equation of motion 

F)V 

V^V^ - — = -kQT (m) . (13) 

For further details see for example 

The FLRW metric in the Einstein frame is 

ds 2 = 5 2 (t?) (-drj 2 + tfydarW) (14) 

where the coordinates are the same in the Einstein frame as in the Jordan frame - this is possible because we are 
employing conformal instead of coordinate time. 



III. THE EQUIVALENCE BETWEEN THE FRAMES 



The interpretation of this "equivalence" has a tangled history in the literature. Usage can be generally separated into 
two camps (see for example who identify a number of works to that date with one camp or the other, and Q which 
sets out clear definitions of "equivalence"): it can be taken to imply that the physics in both frames is identical, or it 
can be taken to imply that a system set up in the Jordan frame can be solved in the Einstein frame as long as it is 
transformed back to the Jordan frame for interpretation. The former case relies on us clearly stating what "physical 
equivalence" means; in [?} (and, similarly (59[) the authors take the reasonable definition that the observables should 
remain the same, so long as the correct length and time-scales are employed. The latter case requires us to define the 
concept of the "physical frame" , the one in which observations should be taken. This would be the frame in which 
it is natural to define our theory; if we are motivated, as in f(R) gravity, by a modification of the Einstein-Hilbert 
action then the Jordan frame would be the physical frame. If instead we were motivated, as in coupled quintessence, 
by exotic couplings between a scalar field and the matter sector, then the Einstein frame would be the physical frame. 

Particularly clear discussions of this issue are found in [l|, Q and for some time it seems generally agreed that 
the "equivalence" is mathematical in nature, but since it occasionally reappears in the literature (see for instance 
0, H, l59 - l68j for some examples since 2004) we briefly re-address the question here. 

"Physical equivalence" between the frames implies that the general behaviour of solutions in the two frames would 
coincide. While the causal structure is naturally preserved by the conformal transformation, it is straightforward to 
find models with very different behaviours. For example, in polynomial gravity with f(R) = aR n one can find [69( a 
dust-dominated cosmology with 

a oc ry-W^-s) ; g K ^n-S) / {m-S) (15) 

in the Jordan and Einstein frames respectively. The deceleration parameter 5 is then 

1 /3-2n\ _ 2ra-3 
*=2-(— )' * = — • (16) 

The universe in the Jordan frame is then accelerating if n > 3/2, but is accelerating in the Einstein frame if n G (3/2,3). 
Choosing n = 4 gives q = —5/8 and q = 5. In the first case the universe is accelerating, while in the second it is 
decelerating. In this strict respect, the frames are clearly not physically equivalent. This example is similar to that 



4 However, it should be noted that more recent coupled quintessence models such as that in |57|| tend to couple the scalar field only to 
particular species of matter, and do not easily admit a Jordan frame representation. 

5 This is defined with respect to the coordinate time t as q = — a(d? a/ dt 2 ) / a? . In conformal time this becomes q = —'H/H 2 . As such, q is 
not the deceleration parameter q transformed into the Einstein frame, but is the deceleration parameter one finds assuming the Einstein 
frame to be "physical". The discrepancy arises because, unlike the conformal time, coordinate time is transformed by the conformal 
mapping. 
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in [67| ; extreme forms are presented in |62| . |63| where the authors find a phantom behaviour in one frame but a 
non-phantom frame in the other, and in [66] where a model based on a non-minimally coupled scalar field exhibits a 
bouncing behaviour in the Jordan frame but not in the Einstein frame. An interesting example from outside cosmology 
is presented in , which studies Birkhoff 's theorem and argues that the Einstein frame representation is stable while 
the Jordan frame is not. 

In the previous section the two frames are shown to be mathematically equivalent at the level of the action. It is 
quick to confirm that the field equations in the two frames are also equivalent to one another, in that a solution in 
the Jordan frame maps directly onto a solution in the Einstein frame and vice- versa. Consider the Einstein equations 
in the Einstein frame, 

- \g^R = k 2 (f^ + f w) . (17) 

Employing the transformations (j6)) and the form of the scalar field stress-energy tensor (| 12[) , it is straightforward to 
convert this to a Jordan frame equivalent, 

R a p ~ \g a pR - 2V Q V (3 w + 2g aP V»V^ - 4V a wV^ + ^^VwV^ = k 2 0- 2 T Q(3 + n 2 V(0)n 2 g a p. (18) 

Using now that F = O 2 > and 2il' i K 2 V = il 2 R — f(R) produces the field equations in the Jordan frame ([2]). Since 
the transformation is invertible, the converse naturally follows. Therefore 

G^, v = k T^ v s — v G^u — k T^ v (19) 

so long as F > is always obeyed. This implies that any solution of the field equations in the Jordan frame induces 
a solution in the Einstein frame, and vice-versa: the two frames are equivalent, even if only one frame is "physical" 
in the sense defined above. See for example [H, 0, EH, [H[ f° r other discussions and interpretations. 

For our purposes it suffices to take the minimal view that the equivalence is mathematical and that observations 
should be made in the Jordan frame. This satisfies the interpretations of both camps; at worst, we have performed 
an unnecessary transformation before measuring the observables. 

IV. PERTURBATIONS AND GAUGE ISSUES IN f(R) COSMOLOGY 

A. Transforming Perturbations between Frames: Gauge Ambiguities 

Applied to a linearly perturbed metric and writing F for the background quantity, the conformal transformation 
between the frames is 

9m» = 9$ + V = ( F + SF ) (< 5 + V) • (20) 
Assuming a flat FLRW background for illustration, all components of the metric imply 

a = ^fFa^U = U + )-^-. (21) 
z b 

The metric perturbation transforms as 

V = ^ V + 9$&F- (22) 
The perturbation can be expanded into scalar, vector and tensor modes [3t| by 

h^dx^dx" = a 2 {r)) (-2$d?7 2 + (d t B - S*,) drjdx 1 + f-2*<% + 2d i d j E + 2d {t F 3) + h\J^ dx l dx j ^ (23) 
with d l Si = d l Fj = d' l h\J = S^h^ = 0. Applying the transformation leads ultimately to 

with all other metric quantities invariant. The transformation then induces both a lapse and a scalar 3-curvature 
regardless of the gauge chosen for h^ v . If for instance we employ uniform curvature gauge in the Jordan frame 
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= E = 0), giving us two scalar degrees of freedom $ and B, then upon transformation into the Einstein frame we 
apparently have three scalar degrees of freedom: <!>, B and 'J. Obviously these are not genuine degrees of freedom, 
being linear combinations of the two Jordan frame degrees of freedom. However, if the field equations are taken at 
face value, one should treat them as three. Erroneously employing the uniform curvature gauge field equations in 
the Einstein frame would in principle induce a gauge mode that would contaminate the results in the Einstein frame. 
Likewise, if we employed synchronous gauge in the Jordan frame, with $ = B = and the degrees of freedom ^ 
and E, upon transformation into the Einstein frame we would have three apparent degrees of freedom: VP, E and the 
redundant lapse $. Clearly if care is not taken the gauge is being tangled in the transformation between Jordan and 
Einstein frames - only gauges possessing both a lapse and a curvature, such as Newtonian gauge, are unaffected. 

The transformations of the matter and scalar field between frames provide a further constraint. Let u^u v Tj^f^ = 
PmO- + <5m) and u^ M ^ = (l/a)(l — 9*i»(Af))> and expand the scalar field ip = K<p as ip — > ip + Sip. Then 

p M = s M = S M - 2—, ip = — \r\nF, 5tp = — -— (25) 

with the velocity remaining invariant. The density perturbation also changes, implying that if we chose uniform density 
gauge in the Jordan frame we would no longer be in uniform density gauge following the conformal transformation. 

Note that this issue only affects scalar perturbations and not vectors or tensors. It is also important to note that 
this is not unique to cosmological models; any system which is separated into a background metric plus perturbations 
will be at risk. 

The cause of this is that the group associated with the conformal transformation is more tightly constrained 
than that of GR. The gauges that lie within this group are those with $ ^ 0, ^> ^ and Sp ^ 0; synchronous, 
uniform curvature ("flat") and uniform density gauges do not belong to this group, while Newtonian ("longitudinal" 
or "Poisson" ) and forms of comoving gauges do. Practically, there are three straightforward approaches to dealing 
with the problem: 

• One could choose to work directly with the additional redundant degree(s) of freedom in the Einstein frame, 
keeping those terms in the field equations. 

• One could "refix" the gauge after a transformation. Transforming uniform curvature gauge between the frames 
induces a lapse; applying a gauge transform into Newtonian or back into uniform curvature gauge will reabsorb 
these, explicitly leaving a system with two degrees of freedom. 

• More physically, one could work in a gauge lying within the restricted group. The most obvious example is 
Newtonian gauge, while other options would include a subset of the comoving gauges with non-vanishing lapse 
and curvature. 

It must be emphasised here that this is only an issue when the conformal transformation is applied - while in 
either frame, one may choose to work with any gauge, as normal. For instance, one should be able to transform in 
Newtonian gauge into the Einstein frame, and then convert to uniform curvature gauge to undertake a calculation, 
before transforming back to Newtonian gauge and transferring the result into the Jordan frame. As a result studies 
in Jordan frame that employ a gauge outside of the restricted class (such as 32] ) remain entirely valid. 

While employing a system with apparent redundant degrees of freedom might seem a bit odd, it can be easily 
justified. By the equivalence of the field equations (|19l) any solution valid in the Jordan frame is equally valid in the 
Einstein frame. A corollary is that, as long as one employs the correct Einstein equations in the Einstein frame, a 
Jordan frame solution will always be valid. However, these correct Einstein equations are not necessarily those that 
would be expected, and must additionally include the additional terms induced by the transformation. For example, 



the gauge-unfixed Hamiltonian constraint 36] in the Einstein frame is 



3H(^> + H^ -d l dt (V + H (e-B^ = -^n 2 a 2 6p. (26) 

If we employ uniform curvature gauge in the Jordan frame, with VP = E = 0, then after transforming into the Einstein 
frame we would naively evolve the system using 

37i 2 $ + Hd l d t B = -^K 2 a 2 Sp. (27) 

However, while we still have E = 0, equation ([24)) implies that 1/0. The correct version of the Hamiltonian 
constraint is then 

3H + - d i d l (* - UB} = -^K 2 a 2 Sp. (28) 
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The validity of this system is guaranteed by equation (fT9|) . 

This might be unpalatable on practical grounds: a great advantage of working in the Einstein frame is to leverage 
existing codes which are written in a set gauge, and evolving the extra degree of freedom will involve additional coding 
and testing. One does, however, remain free to apply a gauge transformation into a more standard gauge. It is this 
which we refer to as "refixing" . The validity of the redundant system in the Einstein frame equally guarantees the 
validity of the refixed system given a well-defined gauge transformation. As such, gauges such as uniform curvature 
and uniform density gauges can be seen as being clean when transforming between the Jordan and Einstein frames - 
and indeed are, so long as either the redundant system is evolved, or is reabsorbed following a gauge transformation. 
Failure to take one of these two steps will result in errors. 

At this point it should again be emphasised that CMB codes are commonly written in synchronous gauge [H], HI, 
EOj EJ ■ The transformation into synchronous gauge contains an arbitrary function of space |36j . While in conformal 
synchronous gauge the ambiguity can be removed by locking synchronous gauge to the velocity of cold dark matter, 
this additional step must be taken if one is intending to use such a code in the Einstein frame. If one takes the 
Jordan frame to be the physical frame (in which the Hubble rate, matter densities and magnitude of perturbations 
are set), then the initial conditions must be set up in the Jordan frame and then transferred into the Einstein frame. 
To employ these in a Boltzmann code a further transformation into synchronous gauge is required regardless of the 
gauge chosen in the Jordan frame, and must be treated with the appropriate care. 6 

The final option is to work in a gauge lying in the restricted group. This group requires the existence of a lapse, a 
spatial curvature and matter perturbations. In the Jordan frame we have four metric and two matter scalar degrees 
of freedom, $, B, ^, E, Sp and v, of which we can remove two with gauge transformations. Of these, $, ^ and Sp are 
non-zero in a gauge belonging to the restricted group, leaving us with the freedom to remove B, E and v. In vacuo, 
this limits us to Newtonian gauge alone, while for systems with matter sources we have greater freedom and can also 
select comoving gauges with v = and either E ^ or B ^ 0. 

In the remainder of this paper we illustrate this issue with a simple example system. 



B. Gauge Issues in Vacuum Polynomial Gravity 

In this section we focus on polynomial gravity, with 

f(R) = aR n . (29) 

The background cosmology of this theory has been considered in detail in (69[. These authors found the exact 
background solution both in vacuum and in the presence of matter fields, 

a K t (l-n)(2„-l)/(„-2)_ (30) 

The exponent diverges when n — > 2, suggesting that it is a particularly interesting case. 
The situation is clearer in the Einstein frame. The scalar field potential ([12"]) becomes 

^ ) = 2a l/(n-D „n" fr-l) ( 2 ^5l ^ ( 31 ) 

which is constant for n = 2. For particular initial data or at late times we would then expect (dtfi/dt) 2 <C V(<f>), and 
the spacetime in the Einstein frame will become de Sitter. For a field with ip = 0, equation (|2Tj) implies that the scale 
factor retains its character - an exponential expansion in the Einstein frame will map onto an exponential expansion 
in the Jordan frame. When n = 2, the solution in the Jordan frame is therefore de Sitter. This has been noted 
before in [7(3], where n = 2 is shown to correspond to w c s = — 1, and in 69] where it is stated that a oc (n — 2) = 
corresponds to de Sitter space. 

We take n — 2 and for simplicity we work on superhorizon scales in vacuum and employ conformal time to ensure 
that our time coordinate remains invariant under the conformal transformation. This simplifies the treatment of the 
perturbations considerably. 



An interesting point that can be made here is that the conformal transformation does not necessarily preserve adiabaticity of initial 
conditions - that is, even if one sets adiabatic initial conditions in the Jordan frame, the introduction of an effective scalar field in the 
Einstein frame may induce a level of isocurvature, which should be taken into account. This issue is left to future study. The authors 
are grateful to David Mota for highlighting this point. 
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Background 



The Einstein field equations in the Jordan frame are 



a\ 2 a 2 n n d „aR „ a ( a\ 2 R a R 



. R-3-+3— - = 0, -+ - — 5 ^ — -rR = 0. (32) 

4 a a-R a \a J R aR 4 

The solution (|30[) is in coordinate time. A solution for this system in conformal time is 

a(r 1 ) = (l + H Q ( Vo -ri))- 1 , H = aH (33) 



where we have normalised such that a(r/ Q ) = 1. Conversely, writing the scalar field and energy density in units of k, 



ip = K<f> and i?^ = K 2 pd>, the Friedmann and continuity equations in the Einstein frame are 



2 



:-!(-) = a 2 R <t .. /?,, + >!«(/?. 

where 



fa + iV) = (34) 



i?0 = ^ 2 + « 2 ^), h = R* - 2« 2 Vm V-(0) = ^-. (35) 



The solution (|3"3"|) transformed into the Einstein frame is 



1 TT 2 



i> = -7^ In (24aif^) , w = -1. (36) 



l + F (»J0-»?) 2Q 

The effective scalar field in the Einstein frame is frozen on a constant potential; as in the Jordan frame, this is de 
Sitter space. 

This solution corresponds to the case where ip = const., which reduces the effective scalar field energy and pressure 
to constants. The system in the Einstein frame admits a more general case, corresponding to a field which is initially 
rolling in the potential before slowing through Hubble friction. The continuity equation for the scalar field can be 
written 

d /_«- \ 3a 5 



<9a V 9 J 4a 
which can be solved to yield 

= £ + & (38) 

If the initial conditions are tuned such that A = 0, or when O 1, then ip = const and this solution is the de Sitter 
limit found above. There is also a "dynamic limit" where the density decays rapidly with time, R^ = A/a e , with 
a = (2y/ A/3(j)o — i])) 1 / 2 , which holds when a -C (SaA) 1 ^ . An implicit general solution for a is 



- t/ '16a 2 \ 1/6 V3tt 2 V24^ fl 17 8Aa , 
V= *{—) (iW ^U'2 ; 6 ; -— ) (39) 



where 2-7 r i(a; b; x) is a regular hypergeometric function. In the limit a oo this tends to de Sitter space. (Note that 
this solution is not normalised such that a(ry ) = 1-) 

For our purposes it is sufficient to work in the de Sitter limit. 



2. Linear Perturbations in the Restricted Group 



In this section we find solutions for the perturbed vacuum spacetime in Newtonian gauge and demonstrate that the 
system is trivially valid in the Einstein frame in both the Newtonian and the uniform curvature gauges. Newtonian 
gauge lies in the restricted class of gauges, while uniform curvature gauge does not. The restricted gauge group 
consists of the usual gauge group in relativity with the additional conditions $ ^ and ^ ^ 0. The uniform density 
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and uniform field gauges are other examples of gauges satisfying these conditions. On large scales where spatial 
derivatives can be neglected the vacuum perturbation equations [9j are 

dF 

F($-*) + —5R = 0, 

6HF ( * + H$ ) + 3— HSR - 3H— [ — ]8R- iH—SR + 3F ( + * J = 
V / oR or) \ oR J oR V / 

where it is understood that F and R are background quantities. Inserting our model with F = 2aR and using 

5R = -% ( 2-$ + H$ + 3-H* + * ) (41) 
a A \ a J 

one can find that 

\J>9 \|/, \Jr 9 \J>, 

tf = tf 1 + _2 + 3 $ = 11-*-^ (42) 

a n .3 



is a general solution. Therefore 



SF ( * 3 . 



This solution is transformed to the Einstein frame using the transformations defined in equations (|24p and (|25l) . 
giving $ = ^ and 

- = ^ M2; ^ = 1 / i+ (24a)3/^ 3 y 



Q 

Then ^ = at a background level, implying 8aR,j, = 1 and = <5P0 = 0. The Newtonian potentials are equal, 
as should be expected since in the Einstein frame, unlike the Jordan frame, the anisotropic stress vanishes. The 
Klein-Gordon equation, Hamiltonian constraint and shear evolution equations [36j j in this frame are 



H ($ + "H*J =0, $ = <5V + 2^5^ = (45) 

which are clearly satisfied. This illustrates the expected equivalence between the frames: a system perturbed in 
Newtonian gauge can be evolved in either the Jordan or the Einstein frames, to yield a consistent result in the Jordan 
frame. The same holds for any gauge that lies within the restricted group. 

A scalar quantity transforms under a gauge transformation ^ — (a, d l /3) as SAp = 5 A + Aa, and so is gauge- 
invariant at linear order if it is constant on the background. As a result, the field perturbation is gauge-invariant 
because "0 = 0, and will therefore always satisfy its equation of motion. The transformation from Newtonian into 
uniform curvature gauge, defined by choosing typ = Ep = [36[, is 

$ F = $ + * + — U , B F = -^, (46) 

where a subscript F denotes uniform curvature gauge. Uniform curvature gauge was used in, for example, [71]. The 
lapse and shift in uniform curvature gauge are then 

24<xff2 

$ F = 0, B F = r T - 2 -1'2 47 

The Hamiltonian constraint and shear evolution become 

3-H 2 $ F = 0, B F + 2HB F = (48) 



which are clearly satisfied. 
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3. Refixing a Gauge 

Consider the system in the uniform curvature gauge in the Jordan frame. Transforming the Jordan frame solutions 
into uniform curvature gauge gives 

\|/, 1 / $,\ 

^ = -^-4-1, B F = + + . (49) 

Since the background Ricci scalar is constant, the perturbed Ricci scalar is gauge-invariant and 5F/F is given by 
equation (|43[) . Transforming into the Einstein frame gives 

$ F , 2 = -3(24a) 3 / 2 iJ 3 |f , * F , 2 = -* x - (24a) 3 / 2 if 3 |^ 

Pur / vi/ * \° (50) 

Bf,2 = ( *i + %/24^o-^ + (24a) 3 / 2 J ff 3 5 ) . 

which evidently differ from (|47[) . In particular, since the lapse is non- vanishing the Hamiltonian constraint is satisfied 
only if ^\ = 0, while the stress evolution equation is only valid if = \&3 = - neither of which are acceptable. 

To illustrate, consider a model in units where a — Hq — 1, and set the initial conditions ^> = 1, ^ = 0, \5 = 1/2 in 
the Jordan frame at a — 1. Then we find 

7 11 

*i = -, ^2 = --, = — ■ (51) 
1 6' 4' 3 12 V ' 

Transferring into the Einstein frame, we find that a = 1 corresponds to a — y24. On this time-slice equation 
(f50|) implies that Bp^ia — \/24) = — 1- Equations (f4"5)l then implies that 

24 

B F ,2 = -^o ( 52 ) 

is the "solution" in the Einstein frame. However, equation (14"T1) shows that 

24 fi 

B F = - — * 2 = — . (53) 
a z a 

The discrepency is of the same order of magnitude as the perturbation itself. Even worse, transferring the lapse gives 
$F ; 2(a = \/24) = —1/4 on the initial hypersurface but equations (14"51) enforce <J>f = at all other times. 

This illustrates the main result of our paper: if you naively transfer a system in an unsafe gauge between the 
Jordan and Einstein frames, you will inevitably introduce errors of the same order of magnitude as the perturbations 
themselves. 

In the previous section we argued two approaches to dealing with this: refixing the gauge, and employing a redundant 
set of field equations. Consider first refixing the gauge to uniform curvature gauge by reabsorbing the spatial curvature 
with the general transformations 

d I 9 \ 9 

$RF,F = $F,2 + *F,2 + 7T" , BrF.F = B F . 2 - -44 (54) 

or] y n J H 

implying ^rf^f = and -Brf,f = — 2AaH^ 2 / <4 in agreement with (|4"T)) . Likewise, refixing to Newtonian gauge, 

i _ . \jr 2 

$RF,N = $F.2 + B F . 2 + UB F<2 = V24aH — , * RF ,N = *F.2 ~ ft -Bp, 2 = $RF,N (55) 

a 

in agreement with (14"4"|) . These illustrate the validity of refixing the gauge to reabsorb the unexpected degree of 
freedom. 

We can also consider the unfixed system in the Einstein frame. We have three scalar degrees of freedom: $p, *f 
and B F . Retaining these three variables, the vacuum field equations on large scales (3(| are 

3ft (* F ,2 + ft^F.a) = 0, B F ,2 + 2ft5 F ,2 + <E>f,2 - *f,2 = 0, * F ,2 + 2ft# F ,2 + ft$F,2 + (ZH + ft 2 ) $F,2 = (56) 
noting that the shear a = E — B. These field equations are also clearly satisfied by equations (|5D|) . 



11 



V. DISCUSSION 

In this paper we have discussed a technical subtlety of the conformal transformation between the Jordan and 
Einstein frames that to the best of our knowledge has not been highlighted before: naively transferring a perturbed 
system between the frames tangles your choice of gauge. It should be emphasised that this effect is calculational, 
not physical, and that if sufficient care has been taken no previous results are changed by this. It is also extremely 
important to note that this issue is not restricted to f(R) gravity or to the study of cosmological perturbations, and 
that it potentially arises whenever a perturbed system is conformally transformed regardless of the background metric 
or the model of gravity. f(R) theory and vacuum cosmology provide a useful, straightforward illustration. While we 
have concretely demonstrated the issue only for this vacuum f(R) system, in principle it arises generically, although 
demonstrating this explicitly is beyond the scope of this paper. 

The properties of the conformal transformation restrict us to a particular set of gauges which we refer to as the 
"restricted group" : gauges that possess a lapse, a spatial curvature, and a density perturbation. This restriction 
only applies when the transformation is performed. At all other times one is free to work in any convenient gauge. 
However, we additionally argued that there are alternatives to working in the restricted group - one can instead 
choose to "refix" the gauge, or to evolve the system using the redundant Einstein field equations. 

We have illustrated our arguments with a simple vacuum system, showing that there are no ambiguities introduced 
by equations (|24l |25|) , so long as one takes such care. If this is not done and the results of a transformation from a gauge 
outside the restricted group are interpreted as if they themselves lie in that gauge, the system will be inaccurately 
specified. 

Perhaps the most important place where this issue will arise is in the treatment of initial conditions. For instance, 
Boltzmann codes are frequently [l4l . l38l . l40l l4l| written in synchronous gauge, which lies outside of the restricted 
group. Interpreting the Jordan frame as physical, the initial conditions must be set in the Jordan frame. However, 
regardless of the gauge chosen, the transformation into the Einstein frame induces a spurious lapse function which 
must be reabsorbed if the results of the calculation are to make any sense. This itself introduces an additional 
problem: the transformation to synchronous gauge is not fully specified, and introduces an arbitrary function of 
the 3-coordinates which must itself be removed with care. The alternative is to include the redundant degrees of 
freedom in the Boltzmann codes, which requires much additional effort. A careless study without attention to these 
issues would either leave the system with the spurious lapse, or risk introducing a gauge mode. Some Boltzmann 
codes [37l HH, H3] include modules written in Newtonian gauge, which lies within the restricted group, and studies of 
extended gravity that employ these codes remain consistent. However, COSMICS is very outdated, while cmbeasy 
and CLASS are programmed in CH — h and C respectively; given the widespread use of Fortran in cosmology, its large 
codebase and the extensive testing it has undergone, CAMB remains the dominant Boltzmann code and the study of 
synchronous gauge is necessary. 7 Similar arguments apply to the initial conditions employed in n-body simulations of 
modified gravity performed in the Einstein frame. 

The other time at which the gauge ambiguity becomes important is at the end of a calculation. Taking the 
Jordan frame to be physical, the results of a simulation in the Einstein frame must be transformed back. Given 
that calculations are frequently performed in synchronous gauge then one must either transform into a safer gauge 
before performing the conformal transformation, or deal with the redundant system in the Jordan frame. Otherwise 
one runs a serious risk of contamination. At this stage it seems easier to refix the gauge in the Jordan frame to a 
fully-specified gauge (such as Newtonian or uniform curvature), or to leave the system unspecified. So long as one 
does so consistently, observables will not be affected, since gauge transformations cannot change an observation. 

It is interesting to note that the lensing potential $ + is uncontaminated by the transformation. As such, if one 
wishes to calculate the weak lensing signal on a constant-time hypersurface, this can be achieved transforming from 
the Jordan frame into the Einstein frame in any gauge and not worrying about refixing the system [l^|. For a static 
(or quasi-static) system, this will be a good approximation. However, carelessly evolving the system without due care 
will in principle introduce errors, so that the frame invariance does not imply such a straightforward estimate of the 
integrated Sachs- Wolfe effect. 

In summary, we recommend that authors exploiting the Einstein frame to study extended theories of gravity work 
in the restricted group. This removes the need for additional gauge transformations to reabsorb the extra apparent 
degrees of freedom, or the need to employ redundant Einstein equations. If one chooses the latter options, care should 
be taken to ensure that the result remains consistent. 



7 CosmoLib |39| is also written in Fortran but in Newtonian gauge; however it has appeared too recently to yet have much of an impact. 
However, it could prove extremely useful for studies of modified gravity. 
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